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Abstract 

The effect of non-sphericity of the quantum dot on the eigenvalues and eigenfunctions has been 
investigated for the case of both the finite and infinite barrier. The ground and excited state 
energies have been calculated for prolate and oblate spheroids as a function of eccentricity of the 
spheroid. The analytic wavefunctions giving the admixture of higher angular momentum states 
have been obtained as a function of eccentricity. 

PACS numbers: 73.22.Dj, 78.67.He 



The theoretical and experimental investigation of quantum dots(QDs) continues to attract 
a lot of attention. Recent experiments have demonstrated that the optical properties of 
QDs can be changed by controlling their shaped"—. Non-spherical QDs are found to give 
rise to radiation which is linearly polarized"^ in contrast to that from a spherical dot. In 
considering the formation of dark and bright excitons in QDs, Efros et alS have shown the 
importance of the non-sphericity of the QDs for the study of the stokes shift. The large 
width of optical transitions in QDs has been attributed to the distribution of both shapes 
and sizes of QDs. The lowering of spatial symmetry in QDs by changing the shape can be 
used to uncover an interesting class of physical laws in different quantum systems^. The 
effect of the departure from sphericity on the eigenfunction and eigenvalues of the dots 
have been studied theoretically by several authors^"—. Efros and Rodinaft have calculated 
the electronic energy levels in ellipsoidal QDs using perturbation theory over spherical dots. 
Cantele et al^~— have studied the energy levels of particles confined to an ellipsoidal quantum 
dot, solving the Schrodinger equation numerically and also by variational methods. They 
have studied the optical transition matrix elements and the oscillator strength in different 
size and shape of QDs. Fonoberov et al^^ have calculated electron-hole energy spectrum and 
Coulomb potential energy in tetrahedral QD's using the finite difference method. In this 
work we calculate the energy levels of particles confined in a ellipsoidal QDs using spheroidal 
harmonics. The eigenvalues and eigenf unctions are obtained as expansions in terms of a 
parameter related to the eccentricity of the ellipse generating the ellipsoid. One obtains 
analytic wave functions revealing the extent of mixing of higher values of orbital angular 
momentum as the departure from sphericity increases. This method is readily extended to 
QDs where the barrier at the surface of the ellipsoid is not infinite. This is of great practical 
importance since the QDs are usually either grown on substrates using different epitaxial 
methods or embedded in a matrix of large band gap material such as glass, polymers, liquids 
etc. In all these cases, the electrons reaching the surface will encounter a finite not an infinite 
barrier. The optical properties of QDs are expected to depend on the height of the barrier 
which will be taken as a parameter. 

Eigenvalues and Eigenfunctions of a Particle in an Ellipsoid: A suitable co-ordinate 
system to investigate the motion of particles confined to the interior of a prolate (oblates) 
ellipsoid is the prolate (oblate) spheroidal co-ordinates. We shall first consider the prolate 
case, and the oblate case later. The prolate spheroidal co-ordinates ^, rj, 9 are related to the 



Cartesian co-ordinates x,y,z by the relations: 
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with — 1<?7<1, 1^'C^oo, O<0<27r. In the prolate spheroidal system 
the surface ^ =constant > 1 is an ellipsoid with major axis of length /^ and minor axis of 
length /(^^ — 1)^/^. f is the confocal distance and ^ = -, where e is the eccentricity of the 
ellipsoid. The surface \ri\ =constant < 1 is a hyperboloid of two sheets with an asymptotic 
cone whose generating line passes through the origin and is inclined at the angle d = cos~^ r) 
to the z— axis. The surface ip =constant is a plane through the 2;— axis forming the angle cp 
with the X, ;z— plane. 

The ellipsoid is assumed to be generated by the rotation around z— axis of an ellipse 
having either the major or the minor axis along z. The semi-axis along z will be denoted by 
'a' and the transverse axis as 'b'. Hence the parameter x = (^/b is a measure of departure 
from sphericity, x > 1 is for a prolate ellipsoid and x < 1 is for oblate ellipsod. The 
Schrodinger equation of the particle inside the ellipsoid 
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where m* is the effective mass, The solution of EqnO can be given by ipim = 
Sim{c, ri)Rim{c, ^)e*'"^ which gives the following ordinary differential equations for Sim{c, rj) 

and Rim{c,0 

.2 



and 



d 
df] 

d 
d^ 






+ 



x(c) 



.W 



2 2 
C 1] — 



m 



l-r/2 






m 



Slm{c,1]) = 



Rlm{c,O = 



(3) 



(4) 



where c = ^kf, k = h\^2rn*. In order that ipi^ be single valued m must be an integer. For 

c = the Eqn.Q) reduces to the associated Legendre equation. Therefore it follows that 

Ajm(O) =/(/ + !), I >m and SimiS^, rf) = Pl^irf) where Pl^irf) are the associated Legendre 

polynomials. Sim{c, rj) for c 7^ can be expanded in an infinite series of the fornJ^ 

00 ' 
Sirn{c,v)= E 4"^(c)Pr+.(^) (5) 

r=0,l 



The prime over the summation indicates that the summation is over only even values of r 
when (/ — m) is even and only odd values of r when (/ — m) is odd. The requirement that 
the wavefunction is finite at r^ = ±1 confines the r] - dependence of the wavefunction to that 
of the angle functions of the first kind. Substitution of (0) in Eqn.Q leads to recurrence 
relations involving the coefficients d}'^2{^)-> d^^{c) and (ij.™2('^)- -^^^ small values of c^ both 
Amn(c) and df^^^^/d\^^, where r = n — m + i', can be expanded in power series of the form: 

,ilm 

>^im{c) = E qitc'', -^^ = E P2kC''' (6) 

A:=0 '-™ fc=l 

The procedure to calculate the coefficients q^^ and p2k is outlined in Appendix. 

The radial functions Rim{c,C,) are obtained from the functions Sim{c,^) using the 

relationi^ 

b 

Rac,0 =Je'''^He - 1)"/'(1 - vT^'SiUc,v)dv, (7) 
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These are the radial functions of the first kind which are appropriate for the particles confined 
to spheroid problem and are given by 

Rlm{c,0 = ^TT;^ 2^ d, {C)t ■ Jm+r{(i) (8) 

where pim is a normalization factor and jm+r{cC,) are spherical Bessel function. 

The co-ordinate ^ is related to the eccentricity e of the ellipse of the revolution by ^ = 1/e. 
The radial functions Rim{c, niust satisfy the boundary condition that at the surface of the 
ellipsoid defined by ^ = ^i = 1/e, the wave function vanishes. This leads to the equation 

£ dl"^{c)f ^^'^^'^- j^^Mi) = (9) 

r=0,l '"• 

With the known expressions of d''^{c), Eqn.ijH)) can be solved to obtain c. Since c = |/c/, 
the energy eigenvalues are given by 

2m* 2m* \ f J 
the subscripts / and m of c indicates that the solutions will depend on / and m. The subscript 
n implies that, in general, there will be more than one solution corresponding to different 
zeros of EqnlHl 

The corresponding equations for the oblate spheroid can be obtained from those of the 
prolate spheriod by making the transformation ^ —>■ +iC,,c — > —io^ in Eqn|21and|3] The 



eigenvalues of energy of a particle in an oblate spheriod are determined by Eqn^J where c 
is now a solution of 
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^ di^i-tc)f ^^'^^''^- j^+M,) = (11) 

r=0,l ^• 

The energy eigenvalues EquHUlfor the prolate and oblate spheriod case can be put in the 
form 

^'^''»= (^^^^I^j^^'^W (12) 

where a = a in case of prolate and a = 6 for oblate. e„im(x) are the energies in units of 
{Ti^ /2m*a^) which depend on the ellipsoid aspect ratio x = ^-/^ but not on a or 6 separately. 
The energy levels as functions of x fo^^ both finite and infinite barriers are shown in Fig.l. 
For prolate case (x > 1) the major axis along z is increased keeping the minor axis in the 
transverse direction fixed at the spherical value. In oblate case (x < 1) the minor axis is 
along z and is held fixed at the spherical value while the major axis in the transverse direction 
is increased. The volume of the spheroid generated by rotation of the ellipse around the 
2;— axis, Ivrafe^, can be written as |7r i?^(a/i?) for the prolate and ^nR^ib/ RY for the oblate 
case where R is the radius of the sphere corresponding to the case a = b = R. (We follow 
the convention that the semi-axis of the ellipsoid along z, major or minor, is denoted by 'a'). 

It is observed that the degeneracy w.r.t. m present in the spherical case {x = 1) is 
removed as one moves away from sphericity. In the prolate case as x increases and in the 
oblate case as x decreases the energies decrease which can be easily understood as caused by 
the deconfining effect due to increase in volume. For a particular / and m, energy increases 
as the barrier height increases, the energy for the infinite barrier being the highest. This is 
again due to the fact that for infinite barrier height the particles are confined to a smaller 
volume. The smaller the height of the barrier at the surface, the bigger is the probability of 
finding the particle outside in a bigger volume and hence the smaller is the energy. 

It is observed that the energy of a state for a prolate with a given value of (|) is higher 
compared to that for an oblate spheroid with the same value for (-). This is easily explained 
by the fact that the volume for the oblate case is greater than that for the prolate with the 
corresponding value of (|). For a given x the {l,m) = (1,0) state has higher energy than 
the (1,1) state for oblate spheroid. The situation is just the reverse in the prolate case. For 
the (1,1) state the motion is mostly confined to the xy plane whereas for the (1,0) state 
the motion is along the z— direction. Since in the oblate case the dimensions along the xy 



direction are greater than that in the 2;— direction, it is expected that the (1,1) energy will 
be less than the (1,0) energy in the oblate case due to deconfinig effects. For similar reasons 
the (1,0) state is lower than the (1,1) state for the prolate spheriod. 

In Fig. 2 the difference between the energies of the states for the spherical and spheroidal 
QDs^Esphere — Espheroid) IS plotted as a function of eccentricity e for the ground state of the 
prolate spheriod in CdSe and GaN materials. Keeping the axis in the transverse direction 
fixed at 16A, the axis 'a' along the 2;— axis is increased thereby increasing the eccentricity. 
The increase of the quantity {E sphere — Espheroid) with e implies that the energies of the states 
decrease with increasing eccentricity. This is equivalent to the observation earlier in Fig.l 
that the energies decrease with x- The inset of Fig. 2 shows the variation of the energies as 
a function of e when the volume is held fixed. In this case Espheriod > Egphere and increases 
with e at a faster rate. As the dimensions along the z— direction is increased (giving rise to 
a decrease in kinetic energy), the dimensions along the transverse directions must decrease 
giving rise to an increase in kinetic energy. The net effect is an increase in energy. The 
corresponding results for oblate spheriod case are shown in Fig. 3 and the inset for constant 
volume. Keeping the axis along z fixed at 16A as the transverse axes are increased, the 
energies decrease with e. In this case Egpheriod < Egphere and results are just the opposite 
when the volume is kept fixed. 

In Tables:l and 2 are given the coefficients of jo(c^) and J2{cC,) in the radial wavefunction 
EqnlHlof the ground state (l,m) = (0,0) as a function of eccentricity in the oblate and prolate 
case respectively, when the volume is held fixed. For both the infinite and finite barrier, 
the admixture oi J2{c^) increases as the eccentricity increases with a corresponding decrease 
in the jo{c^) component, but for the finite barrier case the rate of increase is much less as 
compared to that when the barrier is infinite. 

In conclusion, the energy eigenvalues for both the prolate and oblate ellipsoids have been 
calculated as a function of eccentricity for both the finite and infinite barrier case. The 
analytic expressions for the wavefunctions have been obtained. It has been shown that for 
both the infinite and finite barrier, the probability of states with higher values of angular 
momentum increases as the eccentricity increases, the rate of increase being smaller for the 
finite case. In the context of recent experimental findings^"— , the results presented here 
illuminate the important role of shape in addition to the size of QDs on the possibility of 



tailoring their electronic and optical properties. 
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Table I. 



Eccentricity 
(e) 


Infinite barrier 


Finite barrier 


coefficient of 


coefficient of 


Jo{cO 


J2(c0 


JoicO 


J2{c0 


.1 


.999 


1.22 X 10-'' 


.999 


1.59 X 10-'^ 


.3 


.989 


1.09 X 10-2 


.998 


1.30 X 10-3 


.5 


.903 


.097 


.988 


1.15 X 10-2 


.6 


.788 


.211 


.974 


2.55 X 10-2 


.7 


.610 


.384 


.942 


5.78 X 10-2 



Table II. 



Eccentricity 
(e) 


Infinite barrier 


Finite barrier 


coefficient of 


coefficient of 


Jo{cO 


J2(c0 


JoicO 


j2(ce) 


.1 


.999 


1.21 X 10-4 


.999 


1.59 X 10-5 


.3 


.985 


1.03 X 10-2 


.998 


1.31 X 10-3 


.5 


.917 


.083 


.987 


1.22 X 10-2 


.6 


.828 


.171 


.974 


2.56 X 10-2 


.7 


.706 


.289 


.946 


5.36 X 10-2 



Table Captions 

Table I : The coefficients of jo(c^) and j2(c^) in the radial wavefunction Eqn^oi the ground 
state as a function of eccentricity for the oblate case. 

Table II : The coefficients of jo(c^) and J2{cC,) in the radial wavefunction i?g|Slof the ground 
state as a function of eccentricity for the prolate case. 



Figure Captions 

Fig.l: Eigenvalues as a function of x = | where a is the intercept along the z-axis and b 
along the transverse axis. a=a and b for the prolate and oblate, respectively. 

Fig. 2: The difference between the eigenvalues for the spherical and spheroidal QDs as a 
function of eccentricity for the ground state of prolate spheroid in CdSe and GaN materials. 
The minor axis is fixed and major axis is increased thereby increasing the eccentricity. The 
inset gives the same for the constant volume. 

Fig. 3: The difference between the eigenvalues for the spherical and spheroidal QDs as a 
function of eccentricity for the ground state of oblate spheroid for CdSe and GaN materials. 
The minor axis is fixed and major axis is increased thereby increasing the eccentricity. The 
inset gives the same for the constant volume. 
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Appendix 

The recurrence relation involving the coefficients d''^(c) obtained by substituting Eqn.(jHl) 
in Eqn.Q can be given by 



{l+m+2+u){l+m+l+u) 2 Jim ( \_i_ 

(,2l+3+2u){2l+5+2u) ^ "'l-m+v+2\'^)^ 



Jim 
"'l-m+y 



(13) 



. il-m+u)(l-m-l+u) 2jlm _ n 

' (2l+2u-3){2l+2u-l)^ "'l-m+i'-2 '-'' 



(r > 0) 



where u = 0, ±2, ±4, Putting z/ = in EqnlT^ one obtains the expansion of \im{c) in 

powers of c^ : 

\ (^\ _ /n I 1\ I 2(l)(l+l)-2m^-l 2 
^ImK^) — H' "I" -^J "I" (2i-l){2/+3) ^ 



(l+rn+2)(l+m+l) 2 '^i-m+2 , {l~m)(l-m~l) 2 '^\'"^^2 
^ (2i+3){2Z+5) ^ d|™„ ^ {2l-Z)(2l-l) ^ d|™„ 



(14) 



In Eqniniboth 



S,''' 



and 



sr 



d'" 



are given by power series in c^ : 



jlm Jim 

<^l-m+2 _ p^^2 ^ ^^^4 ^ ___^^^ ^L^ ^ ^^^2 ^ p^^4 ^ ___ 



rfj- 



Putting z/ = 2 in EqnllHIone obtains 



Jim 
H-m 



(15) 



(/ + 2)(/ + 3)-A.„(c) + ^i^^fS§=|f^c^ 

(;-m+2)(;-m+l) ^2 _ (<+r»+4)(/+m+3) ^2 '^i-m+4 



jlm 
"i-m + 2 



(16) 



(2«+l)(2Z+3) 



(2Z+7)(2«+9) 



di" 



di" 



Substituting Eqn|T3l[ISlin EquUHl noting that '^'^i^'^'^ = 0(0^^) and comparing the coefficients 

l — m 

of c^ and c^ from both sides one obtains p2 and p4. Similarly putting u = —2 in Eqnll3land 
following a similar procedure, pg ^^d P4 can be determined, p^, ps-. and pg, pg.. in EqnlTHl 
are obtained by putting z/ = ±4, ±6 etc. and repeating the procedure. 
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